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ABSTRACT 

Aiming to distinguish two types of progenitors of core collapse supernovae, i.e., one with a core 
composed mainly of oxygen and neon (abbreviated as ONe core) and the other with an iron core 
(or Fe core), we calculated the luminosities and spectra of neutrinos emitted from these cores prior 
to gravitational collapse, taking neutrino oscillation into account. We found that the total energies 
emitted as from the ONe core are < 10^® erg, which is much smaller than ~ erg for Fe cores. 

The average energy, on the other hand, is twice as large for the ONe core as those for the Fe cores. 

The neutrinos produced by the plasmon decays in the ONe core are more numerous than those from 
the electron-positron annihilation in both cores but they have much lower average energies ;< 1 MeV. 
Although it is difficult to detect the pre-supernova neutrinos from the ONe core even if it is located 
within 200 pc from the earth, we expect ^9 — 43 and ^7 — 61 events for Fe cores at KamLAND 
and Super-Kamiokande, respectively, depending on the progenitor mass and neutrino-mass hierarchy. 

These numbers might be increased by an order of magnitude if we envisage next-generation detectors 
such as JUNO. We will hence be able to distinguish the two types of progenitors by the detection or 
non-detection of the pre-supernova neutrinos if they are close enough (< 1 kpc). 

Subject headings: stars:evolution — starsimassive — supernova:general 


1. INTRODUCTION 

Massive stars with the mass Mzams ^ 8 M 0 on the zero age main sequence are supposed to explode as supernovae at 
the end of their lives (Janka 2012; Burrows 2013). The explosion is actually instigated by implosion of the central core, 
which is later inverted and leads eventually to the ejection of outer envelopes and the formation of compact objects 
such as neutron stars and black holes. Exactly how that occurs is still a matter of fierce debates (Janka 2012; Kotake 
et al. 2012 and references therein). The stellar core just prior to collapse will be either consisted of irons (referred to 
as the Fe core in the following) or composed of oxygen and neon mainly (called the ONe core hereafter) (Woosley et 
al. 2002). The initial mass Mzams is the main factor to determine which is obtained in the end: stars on the lightest 
end of the spectrum of massive stars (~ 8 — 10 Mq) will have the ONe cores in the last stage of their lives whereas 
more massive stars will produce the Fe cores (Nomoto & Hashimoto 1988; Umeda et al. 2012). All the Fe cores will 
eventually collapse gravitationally whereas only a fraction of the ONe cores will implode with the rest resulting in 
white dwarfs (Poelarends et al. 2008; Doherty et al. 2015). The masses that separate these different regimes are not 
well determined yet, since these stars commonly experience pulsational instabilities that are accompanied by mass 
ejections toward the end of their lives, which are very difficult to compute numerically (Jones et al. 2013). 

In the quasi static evolutions of massive stars, the carbon burning produces the ONe cores as an ash. If the star 
is massive enough (> 10 Mq), the central temperature becomes high enough to ignite these ash elements to produce 
Si and eventually Fe. In lighter stars, on the other hand, further burnings do not occur because the ONe core can 
be supported by degenerate electrons without any energy generation. The core acquires mass thereafter, however, 
through shell burnings while the outer envelopes shedding their masses through the pulsations and winds (Jones et al. 
2013). If the final core mass, which is determined by these competing processes, is small and the central density does 
not reach the threshold Mcore = 1-367 Mq, the result will be the formation of the ONe white dwarf as an end point of 
the evolution (Takahashi et al. 2013). In the opposite case, i.e., if the core becomes massive enough so that the central 
density could exceed the threshold of 10® ®® g/cm®, then electron captures on Mg would commence (Takahashi et al. 
2013). Once this happens, the core contracts further, accompanied by the rise of density and temperature. When the 
density reaches the critical values of 10®° ® g/cm®, the electron captures on Ne, one of the major elements, start and the 
implosion is much more accelerated, eventually leading to supernova explosions called the electron-capture supernovae 
(ECSNe). When the temperature reaches 10® ® K in the gravitational collapse, Ne and O are ignited at the center, 
which is soon followed by Si burnings, and the deflagration starts to propagate outward to convert the ONe core to 
the Fe core. The latter is actually in nuclear statistical equilibrium (NSE) with iron group elements being dominant. 
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The conversion generates heat and the NSE core becomes much hotter than the ONe core before the conversion. 

In the case of more massive stars, the temperatures are higher and electrons are not degenerate until the formation 
of the Fe core. When the central temperature exceeds 10®'^ K in the Fe core, photo-dissociations of nuclei take place, 
which are actually a gradual change in the nuclear abundance in NSE with lighter nuclei becoming more abundant. 
Since this is an endothermic process, the accompanying pressure reduction triggers the collapse of Fe core, which 
eventually leads to supernova explosions, to which we refer as the Fe-core collapse supernovae or FeCCSNe for short 
in this paper. 

These differences in the stellar evolutions to produce the ONe or Fe core are reflected in the marked differences in 
the density profiles of the progenitors just prior to the core collapse: the ONe cores are in general less massive than 
the Fe cores; the density decreases very rapidly outside the core and the envelope is much more tenuous in the ONe 
cores. As a consequence, the supernova dynamics itself may be different between them. As mentioned already, the 
mechanism of CCSNe has been a long-standing unsolved problem for many decades. The initial implosion of the core 
is reversed by the outward sweep of the shock wave produced by the bounce that occurs for the inner part of the core 
when the central density exceeds the nuclear saturation density (~ 3 x 10^^ g/cm®). Unfortunately, the shock wave 
is not strong enough to get through the core completely and is stagnated inside the core. Researchers are exploring 
the way to revive the stalled shock wave and for the moment the so-called neutrino heating scenario is supposed to be 
most promising (Janka 2012). In this scenario, a fraction of neutrinos emitted copiously from a proto neutron star are 
re-absorbed by the matter below the shock wave and the resultant matter-heating reinvigorates the shock wave and 
produces explosion eventually. The mechanism is notoriously inefficient, however, and hydrodynamical instabilities 
of some sorts are believed to be necessary to boost the heating (Janka 2012). Although there is no consensus in the 
supernova society at present on whether the neutrino heating mechanism really works or not (Burrows 2013), we are 
rather certain that it can produce explosion successfully at least for the ONe core even without the boost by the 
hydrodynamical instabilities. As a matter of fact, detailed numerical simulations demonstrated (Kitaura et al. 2006; 
Janka et al. 2007) that shock revival occurs rather early thanks to the compact core and the tenuous envelope with the 
neutrino signals that lack the sign of the instabilities. The resultant explosion is believed to be weak with the energy 
being ^ 10®® erg. This may not be a problem for ECSNe, however, since SN1054, which produced the Crab pulsar, is 
supposed to be such an event (Nomoto et al. 1982; Tominaga et al. 2013). 

Neutrinos are crucial not only in the mechanism of CCSNe but also in observations as vindicated in SN1987A (Arnett 
et al. 1989). Indeed CCSNe are one of the most important targets of the nascent neutrino astrophysics (Raffelt 2012). 
It is true that the PNS cooling phase following the shock revival makes the greatest contributions to neutrino emissions. 
As a matter of fact, this phase lasts for ^ 10 sec and radiates as neutrinos most of the energy of ~ 10®® erg that are 
liberated by the gravitational collapse and stored as the internal energy in the proto neutron star (PNS) (Sato & Suzuki 
1987; Burrows & Lattimer 1988; Fischer et al. 2012). However, neutrino emissions commence much earlier on. It is 
well known in stellar evolution theory that neutrino emissions dominate photon radiations in the stellar cooling after 
C burning. There are five processes responsible for the neutrino emissions: 1. annihilations of electron-positron pairs, 
2. plasmon decays, 3. photo-pair processes, 4. bremsstrahlungs by electrons and positrons accelerated by nuclei and 
5. electron captures on nuclei and free nucleons. Which process is dominant depends on the density and temperature 
of matter in general (Itoh et al. 1996) but the first two processes are normally the most important in the late phase of 
the massive star evolution. In particular, the pair annihilation is dominant at high temperatures whereas the plasmon 
decay becomes more important at lower temperatures. 

What is important here is that the neutrinos emitted by these processes prior to collapse may be observable if the 
supernova occurs in our vicinity, e.g., within Ikpc. Odrzywolek et al. (2004) calculated the luminosities and spectra of 
the neutrinos emitted by the pair annihilations during the C-, Ne-, 0-, Si-burnings for a 20 Mq progenitor model with 
the Monte Carlo method and estimated the detection events for 6 terrestrial neutrino detectors. They found that the 
mean energies of neutrinos are 0.71, 0.97, 1.1 and 1.8 MeV for the C-, Ne-, O- and Si-burnings, respectively. Assuming 
that the distance to the supernova is 1 kpc, they evaluated the event numbers would be 41 for Super-Kamiokande and 
4 for KamLAND. 

In this paper, we investigate the possibility to distinguish the ECSN from FeCCSN by the observation of the pre¬ 
collapse neutrinos. Employing a latest realistic progenitor model with an ONe core (Takahashi et al. 2013) as well as 
those with Fe cores (private communications with Takahashi), we calculate the temporal evolutions of the luminosities 
and spectra of the neutrinos emitted via the pair annihilations and plasmon decays, the two main emission reactions. 
We also evaluate roughly the expected numbers of detection events both for water Cherenkov detectors and for liquid 
scintillators, taking neutrino oscillations into account. As the representative detectors of these two types in current 
operation, we choose Super-Kamiokande and KamLAND. Although the latter detectors have much smaller fiducial 
volumes than the former, they have better sensitivities at low energies. Note that the average neutrino energies 
predecited by Odrzywolek et al. (2004) are lower than the typical energy thresholds > 5 MeV for these large water 
Cherenkov detectors. In addition to these detectors, we also consider Hyper-Kamiokande, a planned next-generation 
water Cherenkov detector, and JUNO, a larger-scale liquid scintillator under construction. 

The organization of the paper is as follows: the progenitor models for ECSN and FeCCSN are described in Section 
2; our method to calculate the luminosities and spectra of the neutrinos emitted by the pair annihilation and plasmon 
decay is summarized in Section 3 and detailed in Appendix; the results are presented in Section 4, and finally the 
summary and discussions are given in Section 5. 


2. PROGENITOR MODEL 
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We employed realistic progenitor models with 8.4 12 and 15 Mq calculated by Takahashi et al. (Takahashi et al. 

2013). The first one produces the ONe core that will give rise to the ECSN whereas the last two generate Fe cores, 
which will result in the FeCCSNe. The evolution of the 8.4 Mq model is not very different from those of the other two 
until the C-burning. Once a core composed of oxygen and neon is formed, they evolve differently. In the lightest model, 
the central temperature does not reach the value needed to ignite Ne, decreasing initially owing to thermally activated 
neutrino emissions. The ONe core contracts on the Kelvin-Helmholtz time scale as a result. The shell He-burning 
increases the mass of the core, which is now supported by degenerate electrons. As the core becomes fatter, the central 
density and, as a consequence, the chemical potential of electrons also rise there. The 8.4 Mq model is massive enough 
to exceed the density threshold p = 10® ®® g/cm® for the electron capture on ®^Mg. The electron capture enhances the 
contraction of the ONe core. Since the electron captures are exothermic, the temperature increases rather rapidly in 
this phase. When the central density reaches p = 10®° ® g/cm®, ®®Ne starts to capture electrons, further accelerating 
the core contraction. When the central temperature reaches 10® ® K, the 0-|-Ne burning takes place. Then the 
temperature rises quickly and the nuclear statistical equilibrium, or NSE, is established at T > 5 x 10® K. The NSE 
region expands at ^ 10® km/s as the burning of 0-|-Ne propagates outward as a deflagration. Electron captures on 
iron group elements as well as on free protons occur in the NSE region just as in the Fe core of FeCCSNe and finally 
lead to a rapid collapse of the core. 

In the case of more massive progenitors of 12 and 15 Mq, the ONe cores achieve high enough temperatures to 
ignite Ne and subsequently O and Si stably with electrons being non-degenerate, to produce Fe cores finally. The 
masses of these cores grow in time as the Si burning continues and the gravitational collapse is induced normally by 
photo-dissociations of heavy nuclei at T > 10®'^ K. It is a rule of thumb that the more massive the progenitor is, the 
higher temperature they have for a given density. This is particularly so, however, when comparing 12 and 15 Mq 
models with the 8.4 Mq model. In the latter model, the core is much cooler until the 0-|-Ne deflagration produces 
NSE, the fact which has consequences in the neutrino emissions in the pre-supernova stages as shown later. The 
neutrino emissions in the pre-supernova stages from massive stars that produce Fe cores were studied by some authors 
(Odrzywolek et al. 2004). In this paper we are concerned with the neutrino emissions from the ONe cores that will 
yield ECSNe later. We employed the 12 and 15 Mq models mainly for comparison to elucidate qualitative differences 
that the ONe cores may make. 

Figures 1, 2 and 3 present the evolutions of these three progenitor models in different ways. Figure 1 is the 
HR diagram, in which the temporal changes in the luminosities L/Lq and effective temperatures are shown as 
trajectories. It is found that the more massive the star is, the more luminous it is but, otherwise, the evolutions are 
similar to each other. The left panel of Fig. 2 shows the temporal changes in the central densities and temperatures 
as trajectories in the p — T diagram. It is evident that three evolutionary paths are not very different from each other 
initially up to pc 10® g/cm®. As mentioned above, the 8.4 Mq star gets cooled via neutrino emissions efficiently after 
C-burning and the central temperature becomes much lower than those in other two models, for which the central 
temperatures increase continuously up to the onset of collapse. In the lightest model, the temperature rises rapidly 
at Pc ~ 10®® g/cm®, at which the electron captures first on Mg and then on Ne enhance the core contraction and the 
0-|-Ne burning commences at some point, producing the NSE core and heating up the matter in it. After NSE is 
established, the central temperature becomes comparable to or a bit higher than those in the Fe cores of 12 and 15 Mq 
models. The electron captures on heavy elements accelerate the gravitational collapse that leads to core bounce in the 
8.4 Mq model whereas the photo-dissociations of iron groups elements trigger the collapse in the other two models. 
Fig. 3 shows the changes in various quantities at the center as functions of time. Note that the horizontal axis is 
the time to collapse. Since it is not easy to define the times of the onset of collapse unambiguously, we choose them 
simply as the hnal times in the stellar evolution models, at which pc = 10®® ® g/cm® and Tc = 10® ® K for the 8.4 Mq 
model whereas pc = lO® ®" g/cm® and Tc = 10® ® K for the 12 Mq model and pc = 10®® ® g/cm® and Tc = 10®® ® K for 
the 15 Mq model. It is clear again that the 8.4 Mq model is different from the other two qualitatively. This is most 
manifest in the third panel, in which the degeneracy of electrons pe/ksTc is shown. Reflecting the fact that the core is 
much cooler in this model, electrons are strongly degenerate after C-burning. We will find later that this feature will 
cause qualitative differences in the neutrino emissions from these progenitors. It is incidentally pointed out that the 
electron fraction Yc drops rather rapidly in the 8.4 Mq model once the NSE core is formed by the 0-|-Ne deflagration 
and electron captures on heavy elements take place. 

Figures 4, 5 and 6 give the radial profiles of density, temperature, electron degeneracy and electron fraction at five 
different epochs for the three progenitor models, respectively. In the left panels these variables are plotted against 
radius whereas the mass coordinate is used in the right panels. The qualitative differences between the 8.4 Mg model 
and the other two are evident. In fact, we can identify the propagation of the deflagration wave, which sweeps through 
the ONe core from the center, producing the hot NSE core behind. It is also recognized that the electron fraction 
drops only in the NSE region. Figures 5 and 6 show, on the other hand, that the electron capture proceeds in the 
entire Fe cores. In the next section we formulate the neutrino emissivities for the electron-positron annihilation and 
plasmon decay and evaluate them for the profiles given here. 

3. NEUTRINO EMISSION RATES 

^ This corresponds to IO.SMq model in their original paper by Takahashi et al. (2013). After taking into account convective overshooting 
in their stellar evolution calculation, they found that the smaller stellar mass produced essentially the same ONe core at the end of the 
evolution. We have hence decided to adopt this updated number in this paper. 
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Of the five neutrino-emission processes that are supposed to occur in massive stars, i.e. electron-positron pair 
annihilation, plasmon decay, photo-pair process, bremsstrahlung and electron capture, the first two are more important 
than the others in general. Having in mind the applications to stellar evolution calculations, Itoh et al. (1996) obtained 
useful fitting formulae to the energy loss rates for these processes. They also drew a phase diagram in the p — T plane 
to indicate which reaction is dominant for a given combination of density and temperature. According to their results 
and the evolutionary paths of our models shown in Fig. 2, the electron pair annihilation will be dominant for the 12 
and 15 Mq models, which form Fe cores and will hence produce the FeCCSNe whereas the plasmon decay will be the 
most important neutrino-emission process for the 8.4 Mq model, which will result in the ECSN. We will see in the next 
section that this is indeed the case. In order to evaluate the numbers of detection events for the terrestrial neutrino 
detectors, not only the energy loss rate but also the spectra are needed. In the following, we describe the formulae that 
we employed in this work to numerically evaluate the luminosities and spectra of neutrinos for the electron-positron 
pair annihilation and plasmon decay. 


3.1. electron-positron pair annihilation 

When the temperature in the core becomes > 10® K, the number of photons with high enough energy to produce 
electron-positron pairs becomes large. As the temperature increases, these pairs become highly abundant, being in 
chemical equilibrium with photons. Although they annihilate each other to generate photons most of the time, they 
produce from time to time pairs of neutrino and anti-neutrino via weak interaction. 

') i -^e+-|-e“ - >Ve + Ve (1) 

The Feynman diagram corresponding to this process is displayed in Fig. 7. 

The number R of reactions to produce the pair of neutrino and anti-neutrino with four momenta = {E^, q) and 
q'^ = (Ep,q'), respectively, per unit time and unit volume® is given by the following equation in the natural unit 
(c = ft = 1): 

^ &) C) /.- (E.) /,* (K) 64 IMr. (2) 

in which Gf = 1.166364 x 10“^^ MeV“® is the Fermi coupling constant, = {Ee,k) and = {E'^,k') are the 
4-momenta for electron and positron, respectively; /g- and /g+ are the Fermi-Dirac distribution functions of electron 
and positron, respectively; the matrix element squared for this reaction is expressed as 

\Mf = {Cv - CAf (q ■ k) {q' ■ k') + {Cv + Ca? (<? • k') {q' ■ k) + mg2 {Cv - Ca? {q ■ qf . (3) 

In this equation, the coupling constants are given as Cy = 1/2 -|- 2 sin® 6^ and Ca = 1/2 with sin® 9^ = 0.2224 for the 
Weinberg angle 9^. Note that all neutrinos are assumed to be massless, which is well justified for our purposes. 

The expression of R in Eq. (2) can be cast into the following form: 

0^2 

R=^[Pih+fi2h+f^zh] (4) 

(27r) 

In the above expression, /I’s are the following combinations of the coupling constants: /3i = [Cy — C^)^, h = 

{Cy -I- Ca)^ and j3z = Cy — and /’s are the functions of the energies of emitted neutrino E^ and anti-neutrino Ep 
and the angle 9 between their momenta q and q': 

h {E^,E„cos9) = - 27rTF;/A,®(l - cosftf , 

[exp{Ei, -I- Ei^)/T - IJ Ag 

[2ymaxG*l (ymax) (ymin)] [l/max^O (^max) 2/min^O (ymin)] ) 

BT ([Gl (ymax) Gi (ymin)] [ymaxGo (ymax) yminGo (ymin)]) 


+C' [Go (j/max) — Go (ymin)]} > 

(5) 

I 2 = h {Ei^,E^, cos 9 ), 

(6) 

^ 2TvTmlE^Ei,{l-cos9) ^ ^ , m 

^ " [exp(F;,. + E^)/T-l] Ag ’ 

(7) 

Ag = Ejj -\- Eu ‘lEi^Ej} cosO^ 

(8) 

A = Ey -\- Ey - EyEy (3 “h COS , 

(9) 

B = —2EJ^ -I- Ap® (1 -I- cos 9) E^Ec (3 — cos 9 ), 

(10) 


^ R is actually a differential number and Rd^ql2Eud^q'l2Eip is the true number of reactions per time and volume. See Eq. (15). 
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and 7] — (^e E Ep -t- Ep'j /P, tj — /Xg/F: ymax — ^ yrain — Ejj^ij^/T and Gn (y) — Fyi {rj y) Eji {rj y), in 

which the Fermi integral Fn(z) is defined as 


Fr 


r 

{z) = 

Jo 


+ 1 


-dx: 


( 12 ) 


77e denotes the chemical potential of electron; the Boltzmann’s constant is taken to be unity in these and following 
equations. The detail of the derivation of these expressions is given in Appendix. 

The number spectrum for the neutrino or anti-neutrino (denoted by i^i) is expressed as an integral of R over the 
momentum of the partner (referred to as V 2 ) as follows: 


dEp, 


E„ 


{27tY 


d?qp 


(27r)" 2E, 


-R{Ep, Ep,cos9). 


(13) 


Note that ui may be or and the natural unit is employed here. The corresponding energy spectrum is just obtained 
as 

<‘Qe p <*3» 

The total number emissivity is found by further integrating dQ'Y / dEp-^ over Ep-^^ as 


Q'Y = 




d^q d^q' 
{2-Kf 2Ep (27r)^ 2E, 


-R{Ep,Ep,cos9), 


and the corresponding energy emssivity Qe is obtained similarly: 


QY = 


dEp 


d^ 

dEp 


d^q d^q' 
(27r)^ 2Ep (27r)^ 2E, 


■EppR{Ep,Ep,cos9). 


(15) 


(16) 


Finally the energy loss rate for this reaction is given by the sum of the energy emissivities over all neutrino species: 

(17) 


allspecies 

Q= E Qe- 


The left panels of Fig. 8 show the number spectra, Eq. (13), of electron-type anti-neutrino for different combinations 
of pYe and T. In the top panel, the value of pYe is varied around pYf. = 10^° g/cw? with the value of T being fixed to 
T = 10^° K, whereas the latter is changed for a fixed value of the former (pYg = 10^° g/cm^) in the bottom panel. It is 
evident that the emissivity is very sensitive to temperature. In fact, as the temperature increases by ~ 20 %, the peak 
number luminosity becomes greater by an order of magnitude. This is simply because the number of electron-positron 
pairs increases as oc T^. It is also observed that the average energy increases as the temperature rises. The dependence 
on pYe is much less drastic: the emissivity decreases with the value of pYg, since the number of electron-positron pairs 
is reduced in this case. 

For a later comparison, we give here the fitting formula for the energy loss rate Q proposed by (Itoh et al. 1996): 


Qpair = i [{Cv^ + Ca^) + 2 (c(.' + C'a^) + { (Cy2 - CA^) + 2 - C'a) } ^pair 


XgWe ^ /pair. 


(18) 


In which the coupling constants are defined as Cy = 1 — Cy and C'j^ = 1 — Ca, and ypair, g (A) and /pair are expressed, 
respectively, as 


gp,,;, = (10.7480A2 -h 0.3967A°-'^ -h 1.0050) 


- 1.0 


X 1 -h (7.692 X 10^A^ -h 9.715 x 10® A°'®) 

p (A) = 1 - 13.04A2 -h 133.5A^ -h 1534A® -h 918.6A®, 


- 1.0 


-0.3 


/pair 


(oo -h ai^ -h 02 ^^) e 


with 


'P-’- ^3 + 5^y-i+52A-2+63A-3’ 

A= ^ 


5.9302 X 109 K’ 


(19) 

( 20 ) 
( 21 ) 

(22) 
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/ p/Me ^ " ^-1 
gjcrfi j 


(23) 


This fitting formula is supposed to be accurate within 10 % of error in the regime, where the electron-positron pair 
annihilation is dominant over other neutrino-emitting processes. In the left panel of Fig. 9 we compare the energy loss 
rates obtained by the formulae given above (Eqs. (16) and (17)) with those given by the fitting formula, (Eq. (18)) for 
different densities and a fixed temperature (T = 10^° K) and electron fraction (Ye = 0.5). Only the contribution from 
electron-type neutrinos is taken into account in this comparison. It is apparent that they are in excellent agreement 
except at high densities p > 10^® g/cm^, where the electron-positron annihilation is no longer dominant and the fitting 
formula is not accurate. 


3.2. plasmon decay 

Plasmons are quantized collective motions of plasma. They are much like photons as shown in the Feynman diagram 
given in the right panel of Fig. 7, obeying the Bose-Einstein statistics. Unlike ordinary photons in vacuum, however, 
they are massive quasi-particles, having a longitudinal mode in addition to the two transverse modes^. Thanks to this 
property, a plasmon decay to two massless particles is not kinetically forbidden, the fact which is in sharp contrast to 
photons in vacuum. The plasmon decay to a pair of neutrinos. 


(24) 


is one of the main cooling processes in massive star after C-burning. As shown later, it is the dominant neutrino- 
emitting reactions in the ONeMg core until NSE is established. 

We calculate the neutrino emissivity via the plasmon decay, following Braaten & Segel (1993). The detailed deriva¬ 
tions are given in Appendix. The number of reactions, R, to produce the pair of neutrino and anti-neutrino with the 
energies E^, and respectively, per unit time and volume, i?, is given as^ 


^ (1 - cos6») 


V 2 J [1 - exp (Ey-I-£'p)T] 


X \ ^S{fL{E^,E^,cos9)) 


El, -f Ec Ey -\- Ey — Ae 

In —-—-:-h 1 


2A, 


— 2 COS0 {Ey + Ey)^ — 2EyEySm6‘^ + 


Ey Ey Y Ag 


2 {Ey + Eyf 


{Ey + Ey COS 0) {Ey + Ey COS 9) 


3ujp'^ {Ey + Ey) . . ,, 

- ——^- s (/t {Ey,Ey,cos9)) 

{Ey COS9 + Ey) {Ey COS 9 + Ey) 


' EyEy (1 - cos 9) Ey+Ey- Ag 
1 H- „ . . -In- 


{Ey Y Ey) Ag Ey Y Ey Y Ae_ 


1 - 


with the following Jl {Ey, Ey, cos9) and fr {Ey, Ey, cos9)-. 

/l {Ey,Ey, cos 9) = 2EyEy (l - cos6») 

, 2EyEy (1 — cos 0) ( Ey Y Ey Ey Y Ey — Ag 
'- ' -In ^-h 1 


+ 3Wp ^^2 

/t {Ey,Ey, cos 9) = 2EyEy (l - COS 9) 


2A,. 


Ey Y Ey -|- Ag 


3 2 T -^p) 

— —CJ 


2 ^ Ag2 


^ EyEy (1 - cos 9) J Ey -|- Ey - Ag 


{Ey Y Ey) Ag Ey Y Ey Y Ag 


(25) 


(26) 


(27) 


Note that the dispersion relations of the longitudinal and transverse modes are obtained from {Ey,Ey, cos9) = 0 
and /t {Ey, Ey, cos 9) = 0, respectively. 

The number spectrum dQ'^/dEy and total emissivities Q'^ and as well as the energy loss rate Q are calculated in 
the same way as for the electron-positron pair annihilations by Eqs. (13), (15), (16) and (17). The right panels of Fig. 8 
show the number spectra for different pYe (top panel) and temperatures (bottom panel). It is found from the figure 
that the number spectrum is much less sensitive to temperature than for the electron-positron pair annihilation but 
depends more on pYg. It is also evident that the peak energy is considerably smaller in the plasmon decay compared 
with the pair annihilation although the amplitudes are not so different between them. This fact has an important 
implication for observability of the neutrinos emitted by these processes on the terrestrial neutrino detectors. 

The fitting formula to the energy loss rate was provided by Haft et al. (1994) in the following form; 


Qplasma — Y 2Cy ^ Qv ■ 


(28) 


® In some textbooks, only the longitudinal mode is called plasmon. In this paper, we refer also to the transverse modes as plasmon. 
^ See footnote 2. 
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In the above expression Qv is given as 

Qv = 3.00 X (/t + /l) f.y, 

in which the following notations are employed: 

2 1.1095 X 

7 = — 


T2 1 (1.019 X 10-V/Me)" 

fx = 2.4 -I- 0.672 -I- O. 5 I 7 -I- 1.2572 
8.672-h 1.357^ 

“ 225 - 177 -h 72 

X = ^ [17.5 -h logio (2p/Me) - 3 logio T] 
y = ^ [-24.5 -f logio (2p/ye) + 3 logio T] 


fxy — ‘ 


1.05 -h {0.39 - 1.25a: - 0.35 sin (4.5a;) 


(|a;| > 0.7 or y < 0) 


—0.3 exp 


— (4.5a: -I- 0.9)^ } x exp < — 


min (0, y — 1.6 -I- 1.25a;) 


0.57- 0.25a: 


(otherwise) 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


The error of this fitting formula is less than 5 % in the regime, where the plasmon decay dominates in the cooling 
over other processes. We compared the energy loss rate obtained by our formula and that given by Eq. (28) in the 
right panel of Fig. 9, in which the temperature is fixed to T = 10^° K just as in the pair annihilation case. Again the 
contribution from electron-type neutrinos is taken into account in this comparison. It is apparent that they agree with 
each other excellently. Note that the plasmon decay is dominant at pYf. > g/cvp? (Itoh et al. 1996, Fig. 2). 


4. RESULTS 

In the following we present the main results: the number and energy luminosities as well as the spectra for different 
neutrino flavors as functions of the time to collapse. Based on them, we then estimate the expected numbers of 
detection events for different terrestrial neutrino detectors. 


4.1. luminosity and spectrum 

Since stellar cores are not homogeneous, we need to calculate the number and energy emissivities per volume and 
time, (5)v 8'iid as well as the number spectra, dQ'^/dE^, at each point of the core and integrate them with respect 
to radius to obtain the number and energy luminosities, L'^ and together with the observed number spectra, 
dL’^j^ldE,: 


Ln= f QN{r)4:TTr'^dr, 

Jo 

nR 

Le= j QEir)4Trr^dr, 

Jo 


dL 


N 


dE, 


y JO 


dEi, 


(36) 

(37) 

(38) 


Here R denotes the radius of the core surface. We evaluate them at different times so that their time evolutions could 
be obtained. 

Fig. 10 shows the number luminosities of different neutrino species for the three progenitor models as functions of 
the time to collapse. The contributions from the pair annihilation and the plasmon decay and their sum are plotted 
separately for the 8.4 Mq model whereas only those from the pair annihilation are shown for the other two progenitors. 
The left panel displays the results for the electron-type neutrinos and those for the other types are presented in the 
right panel. Since there is no heavy charged lepton in the progenitor cores, there is no difference between the mu-type 
and tau-type neutrinos. To the production of electron-type neutrinos, on the other hand, not only the neutral current 
but the charged current also contribute. The number luminosities of electron-type neutrinos are hence larger by an 
order of magnitude than others in general. This is particularly true of the plasmon decay in the 8.4 Mq model. In fact, 
the number luminosity by the plasmon decay is larger than by the pair annihilation for the electron-type neutrinos 
whereas the former is much smaller than the latter for the mu- and tau-type neutrinos. The difference between the 
neutrino flavors, although not so drastic, are also recognized for the pair annihilation. 

Much more remarkable in the figure is qualitative differences between the ONe core and the Fe cores. Among 
other things, the Fe cores emit a substantial amount of neutrinos from much earlier on than the ONe core. This 













Kato et al. 


is simply due to their higher temperatures as mentioned earlier. The electron-positron pairs are more abundant, 
producing neutrinos copiously. In fact, the ONe core is so cold until the electron capture on Mg is opened that 
electrons are strongly degenerate and, as a consequence, positrons are scarce, making the neutrino luminosity from 
the pair annihilation negligibly small. Unfortunately, the plasmon contribution is also tiny as long as the temperature 
is very low, since only a small amount of plasmons are thermally populated then. It is hence understandable that the 
neutrino luminosity becomes substantial in the ONe core only after the electron captures and the subsequent burnings 
of oxygen and silicon heat it up, which occur less than a second before collapse. Even after the formation of the NSE 
region (T > 5 x 10® K) inside the core, the number luminosity of the ONe core is much smaller than those of the Fe 
cores. As a matter of fact, ~ 10^^ neutrinos are emitted in the 15 Mq model while in the 8.4 Mq model the number 
of emitted neutrinos is only ^ 10®^. Note, however, that the luminosity increases very quickly in the last few hundred 
milliseconds in the ONe core and It is still increasing rapidly at the end of our calculation. At this point the central 
density is 10^^ g/cm® and no more data are available from the stellar evolution calculation. Unlike the Fe cores, the 
ONe core is rather slowly contracting even at this point, with the NSE region being expanding gradually. We hence 
expect the luminosity continues to rise further until neutrinos are trapped inside the core, which will be somewhat 
delayed, since the neutrinos considered here have smaller energies 5 MeV) than those produced by electron captures 
during collapse 10 MeV). 

We show the corresponding energy luminosity evolutions in Fig. 11. Unlike the number luminosity, the energy lumi¬ 
nosity of anti-neutrinos is different from that of their partner neutrinos for the electron-positron pair annihilation. They 
are hence displayed separately in the figure. In the case of the plasmon decay, however, the vector current contribution 
is dominant over other contributions and we ignore the latter completely in this paper. In that approximation, there 
is no distinction between the energy luminosity of neutrinos and that of the partner anti-neutrinos. It is apparent that 
the essential features are nearly the same as those in the number luminosities. The energy luminosities for the Fe cores 
are much larger than those for the ONe core through the entire evolutions up to collapse. Although the luminosities 
are different between the neutrinos and anti-neutrinos in the same flavor, differences are indiscernible in the figure. 
One thing to note here is that the energy luminosity of i/g from the plasmon decay is smaller than the that from the 
pair annihilation, the opposite trend to what we found in the number luminosity (see the left panel of Fig. 10). This 
happens because the average energies of neutrinos are quite different between the two emitting processes, which will 
be demonstrated shortly. Incidentally, both the number and energy luminosities are larger for the 15 Mq model than 
for the 12 Mq model. This is simply because the temperatures are a bit higher for the former than for the latter 
according to the general rule that the more massive the star is, the hotter it is. 

Figure 12 shows the normalized number spectra of different neutrino species, which are evaluated for the pair 
annihilation and plasmon decay individually at different times for each of the three progenitor models. Note that 
upper panels correspond to earlier epochs in this figure. The results for the 8.4 Mq model are plotted in the left panels 
and are limited to the last ~ 200 msec, since the luminosities are too small at earlier times as seen above, it is apparent 
from the figure that the peak neutrino energy for the plasmon decay is ~ 0.1 MeV and is not much changed in this 
phase. These values are substantially lower than those for the pair annihilation. In fact, the neutrinos produced by the 
pair annihilation have ^5—10 MeV on average. This is exactly the reason why the number luminosity of electron-type 
neutrinos from the plasmon decay is larger than that from the pair annihilation whereas the corresponding energy 
luminosities have the opposite order. 

It is also evident that the spectra for the pair annihilation are much broader than those for the plasmon decay and, 
as a consequence, a greater number of high energy neutrinos are emitted in the pair annihilation, the fact that makes 
this process more important from the observational point of view, since the cross section of the inverse /3 decay process 
that is employed for detection by the terrestrial detectors is approximately proportional to the square of neutrino 
energy. We can also recognize that neutrinos have higher peak energies than anti-neutrinos of the same flavor. The 
difference is largest for the electron-type neutrino. The origin of these differences is the disparity between electrons 
and positrons in the cores, with the former being more abundant. 

The middle and right panels in Fig. 12 show the normalized number spectra for the 12 Mq and 15 Mq models, 
respectively. In these plots the plasmon decay is ignored and only those from the pair annihilation are displayed. Note 
that the time spans are quite different in these plots from that for the 8.4 Mq model. This is because the Fe cores 
emit neutrinos longer as mentioned earlier. What is most remarkable here is the fact that the peak energies for these 
Fe cores are much lower than those of the ONe core. As a matter of fact, the peak energies (^1 — 2 MeV) for the 
former is even smaller than the threshold energy Fth Ri 1.8 MeV for the inverse /3 decay. The difference stems from 
the fact that electrons are more degenerate and have larger chemical potentials in the ONe core (see Fig. 3). For the 
same reason, the difference between the neutrinos and anti-neutrinos of the same flavor is much less pronounced in 
the Fe cores although the trend is the same. We will investigate what consequences the features we have seen in this 
subsection may have on the detectability of the neutrinos emitted from the different progenitors in the next subsection. 

4.2. event numbers at detectors 

Based on the results obtained so far, we roughly evaluate the expected number of neutrinos that will be detected in 
the terrestrial detectors such as Super-Kamiokande and KamLAND. Both the water Cherenkov and liquid scintillation 
detectors will observe the inverse-/? decay reaction 


Ve+V 


e+ -I- n. 


(39) 


Pre-supernova neutrino emissions from ONe and Fe cores 


9 


which is induced by the pre-supernova neutrinos when they hit the targets in the detectors. Following Odrzywolek et 
al. (2004), we express the cross section <j{E^) of this interaction as 


a{E,) = ao {gv^ + igA^) E,+p,+ = 0.0952 x 10“^' [cm^], (40) 

in which the emitted positron energy and 3-momentum are denoted by Ef.+ = E,^ — {M^ — Mp) and Pe+ = 

\/ Ee+'^ — rue^, respectively, with = 939.6 MeV, Mp = 938.27 MeV and me = 0.511 MeV being the neutron, 
proton and electron rest masses, respectively; the vector and axial-vector form factors are assumed to be gv = 1, 
gA = 1.27 and the constant uo is given as 


o-Q 


cos^ 0c 

TT 


(1-b A 


« ) 

inner/ 


(41) 


with the radiative correction 
is expressed as 


0.024 and the Cabbibo angle 9c = 0.974. Then the event rate at a detector, r, 


N 

47ri?2 


dE^ 


'-E'th 


(F 1^^^ 


(42) 


in which N and R denote the target number in the detector and the distance between the star and the detector, 
respectively. For simplicity, we assume that the detection efficiency is 100 % above the threshold Flth- The features we 
assume in this paper are summarized in Table 1 for Super-Kamiokande, KamLAND, Hyper-Kamiokande and JUNO. 
For the latter two, the actual numbers will be different® from those listed in the table. Since it is not our intention 
here to give very accurate estimates on the event numbers, we believe that they are good enough. The cumulative 
event number Acum is obtained by the integration over the time: 


A,, 



r dt. 


(43) 


In evaluating Eq. (42), the number spectrum should be modified with an appropriate account of neutrino oscillations. 
We show in Fig. 13 the mixing lengths and the density scale heights along with the density profiles for the three 
progenitor models. We assume here that the neutrino energy is 5 MeV and that the mass hierarchy is inverted. Since 
the envelope of the 8.4 Mq is much more tenuous than those for the other two models and the density gradient is 
much steeper accordingly, the ratio of the scale height to the mixing length becomes smallest for the ONe core model 
at the resonance, which is marked with a star in the figure, just as expected. It is evident, however, that the scale 
height is still larger than the mixing length by a order of magnitude even in that case. We hence assume adiabatic 
oscillations in the following. This is certainly justified for the normal hierarchy, since there is no resonance in the 
anti-neutrino sector. It is noted, however, that the mass and structure of the outer envelope of the 8.4 M© model are 
highly uncertain just prior to collapse. Since the envelope was ignored in the evolution calculation after the completion 
of C burning, we attach it again to the core model, assuming that it is hydrostatic and its temperature and chemical 
composition profiles are unchanged from those at the end of C burning. 

The number spectrum of the electron-type anti-neutrino is then given as follows: 


(dL'^\ 

\ dEjy^ J 


= p 



+ {l-p) 



(44) 


In this expression, the subscript 0 means the original spectra before the neutrino oscillations are taken into account; 
Vx stands for or Vr , both of which we assume to have the same spectra; the so-called survival probability p is given 
in the adiabatic limit as 

_ J lAeip = cos^ 012 cos^ 013 for normal hierarchy, , . 

^ \ I Ae3 P = sin^ 013 for inverted hierarchy. 


with cos^ 012 = 0.692, cos^ 0i3 = 0.9766 for normal hierarchy and sin^ 0i3 = 0.024 for inverted hierarchy (PDG data 
2014). In the above expression, the Uei and Ues are the elements of the unitary matrix U transforming the mass eigen 
states to the flavor eigenstates as 

(46, 


It is normally parametrized by three mixing angles, 0i2, 0i3, 023 and a CP-violating phase, <5, as 


/ Ucl 

Ue2 

C/e3 \ 

/ C 12 C 13 

S 12 C 13 

5136-*^ \ 

U,1 

Up2 

U,3 

1 = ( — S 12 C 23 — Ci2Sl3S23e*'* 

C 12 C 23 — Si2Sl3S23e*'^ 

C 13 S 23 j 

V Url 

Ur2 

Urs) 

\ S 12 S 23 — Ci2Sl3C23e*'* 

-C12523 - Sl25l3C23e'^ 

C 13 C 23 / 


U = 


® For JUNO, as a matter of fact, we just scale KamLAND by a factor of 20. 


(47) 
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with the common shorthand notations of Cy = cos 6ij and Sij = smOij for i,j = 1, 2, 3. 

Fig. 14 show the cumulatice event numbers at Super-Kamiokande and KamLAND as functions of time. The upper 
panel is for the Fe core models whereas the lower one corresponds to the ONe core model. The distance to the star is 
assumed to be 200 pc, the estimated distance to Betelgeuse. In the case of the Fe core models, more than 1 neutrinos 
will be detected a day before collapse both at Super-Kamiokande and KamLAND, which seems to be consistent with 
the preceding paper (Odrzywolek et al. 2004). The total event numbers are also comparable between the two detectors 
although the detector volumes are vastly different. This is mainly due to the sensitivity of KamLAND to low energy 
neutrinos (see Table 1 for the energy threshold of each detector). 

In the case of the ONe core model, on the other hand, the expected event numbers are much smaller than those for 
the Fe core models. This is so in spite of the fact that the peak energy of De is about twice larger for the ONe core 
than for the Fe cores after the formation of NSE core in the former. It turns out that the duration of such neutrinos 
is simply too short in the ONe core. Since the volume of the NSE core is small accordingly, the luminosity itself is 
smaller for the ONe core. 

Both the water Cherenkov detector and liquid scintillation detector may see a major scale-up in the near future. 
For example, Hyper-Kamiokande is a planned next-generation water Cherenkov detector with volume expected to be 
^15 times larger than that of Super-Kamiokande. JUNO, on the other hand, is a liquid scintillation detector with the 
size ^ 20 times as great as that of KamLAND ans is currently under construction. Simply assuming the parameters 
summarized in Table 1, which will certainly different from the actual numbers one way or another, and ignoring the 
detection efficiencies again, we calculated the event numbers for these future detectors. The cumulative event numbers 
are shown in Fig.15 as functions of time. The distance to the source is again assumed to be 200 pc. 

It is found that for the 15 Mq model with an Fe core, the total event number reaches ^ 266 — 864 for JUNO, which 
means that we may observe ~ 10 — 35 events even if the source is located at 1 kpc. Hyper-Kamiokande, on the other 
hand, will observe ^ 28 — 77 events, which is not much different from those for Super-Kamiokande. This is because of 
the somewhat high energy threshold of 8.3 MeV we employ here. If this were as small as the one for Super-Kamiokande, 
the event number would be larger by more than an order of magnitude. The disadvantage of the high energy threshold 
is also vindicated in the expected event numbers for the 12 Mq model. Since the peak energy of neutrino is somewhat 
lower in this model, the event numbers at Hyper-Kamiokande is even smaller than those at Super-Kamiokande. In 
the case of the 8.4 Mq model with an ONe core, the event numbers are expected to be much smaller than unity even 
for these large scale detectors and the distance as small as 200 pc. We are hence forced to conclude that it is highly 
difficult to observe the pre-supernova neutrinos from ONe cores. Put another way, no detection of such neutrinos from 
a nearby CCSN in the future may indicate that it is an ECSN. 

5. SUMMARY AND DISCUSSIONS 

In this paper, we have investigated the pre-supernova neutrino emissions from the ONe-core forming progenitor and 
compared them with those from the more massive Fe-core producing progenitors, aiming to distinguish ECSNe form 
more common FeCCSNe by future observations of a near by supernova. Employing the realistic progenitor models, we 
have calculated the temporal evolutions of the luminosities and spectra of the neutrinos emitted by these stars via the 
electron-positron pair annihilations and plasmon decays. Based on these results, we have then roughly evaluated the 
numbers of detection events for Super-Kamiokande and KamLAND as the representative water Cherenkov and liquid 
scintillation detectors, respectively, in current operation. We have also considered Hyper-Kamiokande and JUNO, 
which may be operational in the near future. 

We have found that the electron-positron pair annihilation dominates over other neutrino-emission processes in 
the Ee-core producing models whereas in the ONe-core generating model the plasmon decay is dominant until the 
NSE core is formed at the center, after which the pair annihilation prevails quickly and the neutrino luminosity rises 
rapidly. Since the contraction proceeds slowly in the latter case, the central density is larger just prior to collapse 
and, as a result, emitted neutrinos have larger peak energies than in the Fe cores, which is advantageous from the 
observational point of view. It is also found that the neutrinos emitted via the plasmon decay have much smaller 
energies O.I MeV) and will be undetected on the terrestrial detectors that employ the inverse decay reaction with 
the energy threshold of « 1.8 MeV. 

Based on these results we have roughly estimated the number of detection events for Super-Kamiokande and kam- 
LAND. We have seen that both will detect tens of neutrinos from the FeCCSNe if they are located at 200 pc from the 
earth whereas neither will be able to detect the ECSNe even at this small distance. If we envisage next-generation 
detectors such as JUNO, these numbers may be increased by an order of magnitude. On the other hand, Hyper- 
Kamiokande may not be advantageous for the detection of such low energy neutrinos unless the energy threshold is as 
low as that of Super-Kamiokande. 

It should be pointed out that the neutrino emissions will continue after the gravitational collapse of the core com¬ 
mences. As mentioned earlier, the core contraction is slow in the 8.4 Mq model even when the central density reaches 
10^^ g/cm^ and it may take another few hundreds msec until core bounce. Since the neutrino luminosity is till rising 
rapidly even at the end of calculation. It is hence expected that the actual event numbers will be increased substan¬ 
tially. Some delay of the neutrino trapping, which will actually terminate the neutrino emissions, owing to the lower 
neutrino energy will be also advantageous in this respect. It is hence important to calculate the neutrino light curve 
up to core bounce, after which ordinary supernova neutrinos will prevail. 

Although we have not considered the plasmon decay for the Ee cores in this paper, it may be comparable to or even 
dominant over the pair annihilation in the late phase according to the right panel of Eig. 2. We hence compare the 
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total emissivities of these two processes just prior to core-collapse for the 15 Mq model in Fig. 16. In the central region 
r < 10^ cm, the plasmon process emits more energy indeed. In the outer part of the core, however, the pair annihilation 
becomes dominant and emit far greater amounts of energy in neutrinos. This is due to the strong dependence of the 
pair annihilation on the temperature (see the lower left panel of Fig. 8). Note also that the plasmon decay is suppressed 
at lower densities (see the upper right panel of Fig. 8). Integrating these emissivities over the radius, we find that 
the total emitted energy from the plasmon process is ~ 10^® erg/s, which is much smaller than that from the pair 
annihilation (^ 10"'^^ erg/s) and justifies our neglect of the plasmon decay in the Fe-core progenitor models. 

So far we have seen that the pre-supernova neutrino emission is qualitatively different between the ONe and Fe 
cores, which fact is reflected in the expected event numbers on the terrestrial detectors. Actual detectability is 
crucially dependent on the background noises at detectors, however. Although it is much beyond the scope of this 
paper to take them into account in detail for each detector and discuss the detection possibility quantitatively, we 
can touch the issue briefly: if we adopt several hundreds events/day as a typical background for Super-Kamiokande, 
neutrinos from FeCCSNe may be detected at 3a if the source is located at 200 pc. The background will be reduced 
substantially if the neutron tagging is enhanced with Gd (Beacom & Vagins 2004). The background at KamLAND 
is low ^ I event/day. The neutrinos from Fe cores may be detected to kpc. For the next-generation detectors, the 
distance may be further extended. In the case of Hyper-Kamiokande, however, the reduction of the energy threshold 
will be more crucial as mentioned above. In this paper, we have considered only two relatively light Fe-core progenitors. 
It is certainly important, though, to study more systematically pre-supernova neutrinos from different progenitors. 
It should be also stressed that the expected event numbers for these models may change by a factor of a few if one 
considers various uncertainties in the stellar evolution calculation. These issues should be investigated more in detail 
and are actually being currently undertaken (Yoshida et al. 2015 in preparation). 


APPENDIX 


In the following we give detailed derivations of the expressions for the reaction rates of the electro-positron pair 
annihilation and the plasmon decay in turn. 

REACTION RATE OF ELECTRON-POSITRON PAIR ANNIHILATION 

The rate of the reaction, R, to produce a pair of neutrino and anti-neutrino via the annihilation of a pair of electron 
and positron, which corresponds to the left Feynman diagram in Fig. 7 to the lowest order, is given by the low-energy 
limit of the Weinberg-Salam theory, which is actually identical to the Fermi’s theory, as follows: 






d^k 


d^k' 


{2Try 2Ee {2TTy 2E, 


- i2ny S^{q + q'-k- k') /,- {E^ /,+ (A') 


x64 


{Cv - Ca? (q ■ k) {q' ■ k') + (Cv + Ca? {q ■ k') {q' ■ k) + {Cl - C\) [q ■ q') 


(Al) 

in which k = {Ee,k), k' = {Ej,k'), q = {E^,q) and q' = {Ey,q') denote the four momenta of electron, positron, 
neutrino and anti-neutrino, respectively. Following Schinder & Shapiro (1982); Mezzacappa et al. (1993), we re-cast 
the above equation into the following form: 


R = 


8G^ 


[Pih + P2I2 + Psh], 


(A2) 


in which the factors are given as Pi = {Cv — Ca)^ , P 2 = {Cv + Ca)^ , Ps = Cy — C\ and the integrals are grouped 
into 


h = 
h = 
h = 


d^k dPk 


Ee E, 
d^k d^k 


-5^ {q + q'-k- k') /,- {Ey /,+ {E’y {q ■ ky 


Ee E, 
dPk dPk 


-5^ {q + q'-k- k') fe- {Ee) fe+ {K) {q ■ k'y, 


Ee E, 


{q + q'-k- k') fe- {Ee) fe+ {K) ml {q ■ q') . 


(A3) 

(A4) 

(A5) 


These integrals are evaluated as follows. We begin with R. Using three of the four (5-functions we can easily 
accomplish the integrals over the positron momentum k' to get 


/ (fly / 9 \ 

-^/e+ {E'e) S'^{q + q' -k - k') = 2/g+ {E^ -b Ap - Ee) 0 {E^ + Ep - Ee) 5 [q + q' - k - mlj , 


(A6) 


in which the Heaviside function is denoted by 0. The remaining integrals over the electron momentum k are performed 
on the spherical coordinates in the momentum space with the volume element written as d^k = \k\EedEed{cos 0e)d(j)e- 
The pe integral is trivial to give a factor of 27r. The integral over 6e can be accomplished with the use of the last 
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^-function. The resultant expression is given as 


h = 


/ 27 r dE, /e- {Ee) fe+ {E, +E^- Ee) " / (1 - cos ef [AEl + Bile + C] , 


in which 6 is the angle between q and q' and is given by 

+ Ei,^ + 2E,^Ep cos 0, 

and A, B and C are defined as 

A = Ep"^ + Ep'^ — EpEp (3 + cos0), 

B = —2Ep'^ + Ep^ (1 + cos0) + EpEp (3 — cos0), 


^ ^17 , 17 a^2 1 77 2 /)2\ 1 /WeAeV 1 + C0S6» 

C = {Ep + EpCosO) --^Ep [l-cos9)--l— —1 


1 — cos 9 

The lower and upper limits of the remaining integral with respect to E^, are given as 


En-iin = max 


i^max = min 


Ep + Ep Ag 

me, ---1 - 


2mp 


EpEp (1 — COS0) 


77 I 77 _ Ep + Ep ^ Ag j ^ 


2mp 


EpEp (1 — COS0) 


Finally, we express the integral over the electron energy Ee with the Fermi integral defined as 


r°° ^n. 

En id) = / —:- rdx. 

Jo 


70 6=^-'' + ! 

Using the following relation for the product of the Fermi-Dirac distributions, 

fe-{Ee)fe+{Ep + Ep-Ee)= ^ 


exp[(£;^ + Ep)/T] - 1 
1 


exp[(£’g - {Ep + Ep) - ^g)/T] + 1 exp[(F;g - fie)/T] + 1 


we can obtain the final expression for Ii as 
27rr£;2£;p2(l-cos6l)^ 


^1 = -^ 


{AT^ {[G2 (j/Tnax) - G2 (j/Tnin)] 
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(AlO) 
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[exp(F;, + Ep)/T) - 1] Ag5 
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with q' = (ye + Ep + Ep) /T, q = yg/T, ymax = E^ax/T, ymin = E^in/T and G„ (y) = F„ {q' - y) - {q - y). 
Similar calculations can be done for the other two integrals to give 


l2 = h {Ep,Ep,cos9 ), 

2TTTmlEpEp (1 — cosy) 


E — — 
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(A13) 

(AM) 

(A15) 


REACTION RATE OF PLASMON DECAY 


The properties of plasmon are derived from the so-called polarization tensor 11'^^', which is calculated field- 
theoretically as 


zn^"(iF) = ^ 


d^k {k ■ K) {Ri^k^ + K‘'k^^) - K^ki^k’' - {k ■ K)^gi^‘' 


2Ee {2 tt) 


{k ■ KY 


ife- (Ee) + fe+ (Ee)) 


(Bl) 


in which K = (w, K) and k = (Ag, k) are the 4-momenta of plasmon and electron, respectively. It is decomposed into 
the transverse (IIt) and longitudinal (IIl) components, which are expressed as 

n'"'^(A") = nT{K)PY''{K) + nL{K)PY’'{K), (B2) 
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where the projection operators and are given as 

Pri^''{K) = 


•0 (/r,z/) = (0,0) 

0 = (0,i) 


(B3) 


\K\ 


Pl^^{K)=1 - 


ojK^ 


= ( 0 , 0 ) 

(/r, z/) = (0, i) . 


■ 


(B4) 


Using the relations P^'^Pl^v = 0; Pt'^Pt^u = 2, P^"Pl^v = 1 between the projection operators, we obtain each 
component of the polarization tensor as 

1, 


nT{K) = -n>^%K)PT^AK) 


(Pk 


-K^ \k 


2 II,|2 


^ {K -kf + 2{k- Kf 


\K 


J 2 Ee { 27 ry 
n.(^)^-|^noo 


{k-KY 


(/e- {Ee) + /e+ {E,)) , 


A f (Pk -K'^EY + 2{k-K) {Ee uj)-{k- Kf 


\K\ 


2Ee {2 ttY 


{k-Kf 


ife- m + /,+ (Ue)) . 


(B5) 


(B6) 


In the relativistic limit, i.e.. Eg ^ \K\ ,ui and rUe —^ 0 and E^ |fc|, which is well justified in the present case, the 
above expressions are reduced to the following: 




3 


IIl{K) = -3uj. 


2e^\K 


,2^P 


^ uj^-\Kf uj ,^+1X1 
1-In ■ 


p 


\K[ 




bJ , bJ ■ 

■ In ■ 


2\K\ bJ-\K\ 
-\K\ 


2\K\ bj-\K\ 


- 1 


in which the plasma frequency ujp is defined as 


bJp^ = ^{kBTf {F,{P+F,{-P). 


(B7) 

(B8) 

(B9) 


The number of the reaction per unit time and volume, R, to produce a pair of neutrino and anti-neutrino via the 
decay of a plasmon is given by the polarization tensor as follows: 


R = 


71 

G 


8 (g/x?!/' + q^qf - Ofi^q ■ q') 2TTCf^K'^ [©(w) (1 -h /b(w)) -H 0(-a;)/B(-a;)] 

X {\iL{K)PY'5 [K^ + e'^ULiK)) + nT{K)PY'^5 {K^ + P^AriK))) 

8K^ 


/?r / h ^ ,,^[^(‘^)-®i-‘^)]{<ipql + q'^q't,-9p<^q-q') 

V2/ [1 - exp (-w/A:bT)J 


3 


X < -- 


2ePK 


2^P 


-3bJr. 


1 - 




2\K 


In- 


\K\ 


-\K\ 


PY''5 {K^ + e^UriK)) 


m 


UJ , UJ ■ 

■ In ■ 


\K\ 


2\K\ bJ-\K\ 


-1 PY''6{K^ + e^ULiK))\ . 


(BIO) 


Note that the dispersion relations for the transverse and longitudinal plasmons are obtained from the i5-functions as 
bJ = ujt{K) and uj = ujb{K), which satisfy the following equations: 


-H e^nT(-ftr) =K^ - - 


3 2 


2-P \Kf 


bj'p — \K''^ lot , + l-^l 

1--m 


lot 


2\K\ ujt-\K\ 


= 0 , 


(Bll) 
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K 


+ e‘^IlL{K) =K^- 3uj/ -^ 


2 r 


\K[ 


WL In^L + IKI 


2\K\ u;l-\K\ 


- 1 


= 0 . 


(B12) 


Employing the conservation law K = — {q + q'), we finally obtain the reaction rate i? as a function of Ei,, Ep,cos9 
as follows: 

^_._Gyi6Cy2 2£;^2Ep2(l-cos6») 


^/2 J [I - ey.^{Ev + Ep)lkBT)\ 

/ Ejj+Ejj E^ + E[, — Af, 

X i-^S{fL {E^,Ec„ cos9)) -In 


Ej/ + Ep + A, 


+ 1 


2 (E -l- E-\^ 

-2 cos 9 {Ep + Epf - 2EpEps\n 0^ + ^ \ {Ep + Ep cos 9) {Ep + Ep cos 9) 


3ujp^ {Ep + Ep)^ 

- —T-^- S (/t {E^ , Ep, cos 9)) 

{Ep cos 9 + Ep) {Ep cos 9 + Ep) 

! o 


' EpEp (1 - cos 6 ») Ep + Ep - Ae 
1 H-tt: - „ ' . -In- 


{Ep + Ep) Ag Ep + Ep + Ag 


with /l {Ep,Ep, cos9) and /t {Ep,Ep, cos9) given as 

Jl {E^,Ep,cos6) = 2E^Ej^ (1 - cosO) 

2 2EjjEy (1 - cos f Ey + Ey Ey + Ey - Ag 


“hScj 


'p 


2Ap Ep + Ep + aE^'^' 


Ag2 

Jt {Ey,Ey,cos6) = 2EyEy (1 - 003 6 *) 

EyEy (1 - COS^) Ey + Ey - A 

1 H- 7 ^- _ . . -In- 


3 2 E Ey) 

“CJn 


■^p 


Ag 


{Ep + Ep) Ag E/j/ + Ep + Ag 


(B13) 

(B14) 

(B15) 
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Figure 1. The trajectories in the HR-diagram of the progenitors employed in this paper. The red, blue and green curves correspond to 
8.4, 12 and 15 Mq models, respectively. 
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Figure 2. The trajectories of the central density pc and temperature Tc for the progenitor employed in this paper. The red, blue and 
green curves correspond to 8.4, 12 and 15 Mq models, respectively. Also shown are the boundaries of domains, in which different processes 
are dominant as indicated (Itoh et al. 1996). 
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Figure 3. The time evolutions of the center central density, temperature, electron degeneracy and electron fraction. The red, blue and 
green curves correspond to 8.4, 12 and 15 Mq models, respectively. 




















Pre-supernova neutrino emissions from ONe and Fe cores 


19 


CO 

'e 

o 

.O) 

Q. 

O) 

o 


hs 

O) 

o 








r [cm] MIMq 


Figure 4. The radial profiles of density, temperature, electron degeneracy and electron fraction in the 8.4 Mq model at different times. 
Different colors denote different times to collapse: t = 0.2338 s (red line), 0.1972 s (blue line), 9.350 x 10“^ s (green line), 6.879 x 10“^ s 
(orange line), 5.982 x 10“^ s (purple line). 
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Figure 5. The same as Fig. 5 but for the 12 Mq model. Different colors show different times to collapse: t = 4.537 x 10^ s (red line), 
9.490 X 10^ s (blue line), 1.344 x 10® s (green line), 6.044 x 10^ s (orange line), 0 s (purple line). 
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Figure 6. The same as Fig. 5 but for the 15 Mq model. Different colors show different times to collapse: t = 8.431 X 10® s (red line), 
1.824 X 10^ s (blue line), 2.099 X 10® s (green line), 9.340 X 10^ s (orange line), 0 s (purple line). 
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Figure 7. Feynman diagrams of the electron-positron pair-annihilation (left panel) and the plasmon decay (right panel). 
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Figure 8. The number spectra of the pair annihilation (left panels) and the plasmon decay (right panels) for different combinations of 
pYe and T. The top panels show the dependence on pYe, in which T is fixed to 10^® K, whereas the lower panels display the dependence 
on T, where pYe is set to 10^® g/cm^. 




Figure 9. Comparisons with the fitting formulae for the pair annihilation (left panel) and plasmon decay (right panel). The temperature 
and electron fraction are fixed to 10^® K and Ye = 0.5, respectively, in both panels. 
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Figure 10. The time evolutions of the number luminosities of electron-type neutrinos (left panel) and of other types (right panel) for 
three progenitor models. 




Figure 11. The time evolutions of the energy luminosities of electron-type neutrinos (left panel) and of other types (right panel) for 
three progenitor models. Top and bottom panels show the results for neutrinos and anti-neutrinos, respectively. 
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Figure 12. Normalized number spectra at different times for 8.4 Mq (left panels), 12 Mq (middle panels) and 15 Mq (right panels). 
Red, blue, green and orange curves correspond, respectively, to Pe and Ve from the pair annihilation and P^/Pr and from the pair 

annihilation. All neutrinos have the identical spectrum after normalization for the plasmon decay as shown with black. For better visibility, 
all the lines but the black one in the left panels are multiplied by the factors indicated. Note that larger t’s correspond to earlier times. 
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Figure 13. Comparisons between density scale heights {Hp) and mixing lengths {Im) for three progenitor models (lower panels). The 
inverted mass hierarchy is assumed and the resonance points are marked by star symbols and vertical lines. The radii of the resonance 
points for the 8.4, 12, 15 Mq models are 1.3 xlO®, 2.95x10^ and 3.53x10^ cm, respectively. In the upper panels we show the density 
profiles for convenience. 
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Figure 14. The cumulative numbers of detection events at Super-Kamiokande (left panels) and KamLAND (right panels) as functions 
of time to collapse. The top panels show the results for the Fe-core models whereas the bottom panels display the evolutions for ONe-core 
model. Note the difference in the scales of the horizontal axis between the upper and lower plots. Red, blue and green lines correspond to 
the 8.4 Mq, 12 Mq and 15 Mq models, respectively and solid, dotted and dashed curves show the results for no neutrino oscillation and 
those for adiabatic oscillations with the normal and inverted mass hierarchies, respectively. 
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Table 1 

The detector parameters assumed in this paper.^ 


Detector 

Mass 

[kt] 

Target number 

N 

Energy threshold 
[MeV] 

Super-K 

32 

2.14x10^^ 

5.3 

KamLAND 

1 

8.47x10^^ 

1.8 

Hyper-K 

540 

3.61x10^^ 

8.3 

JUNO 

20 

1.69x1033 

1.8 


References. — (1) Super-K Calib 2013; (2) KamLAND Calib 2009; (3) Hyper-K Calib 2011; (4) 
JUNO Calib 2014 

^The numbers given here are not very precise and just meant for rough estimate. JUNO is assumed to be a scale-up of KamLAND by a 
factor of 20. We also assume that the energy threshold of Hyper-Kamiokande will be somewhat higher than that of Super-Kamiokande. 
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Figure 15. The same as Fig.14 but for Hyper-Kamiokande (left panels) and JUNO (right panels). 


Table 2 

The expected numbers of detection events for different detectors.*^ 


detector 

8.4 

normal 

Mo 

inverted 

12 

normal 

Mo 

inverted 

15 

normal 

Mo 

inverted 

Super-K 

2.47 xlO"^ 

9.68x10-3 

21 

7 

61 

21 

KamLAND 

1.06 xlO-3 

1.50x10-3 

31 

9 

43 

13 

Hyper-K 

0.30 

0.13 

9 

3 

77 

28 

JUNO 

2.12 xlO-2 

8.03x10-3 

618 

189 

864 

266 


^The source is assumed to be located at 200 pc from the earth. Both the normal and inverted mass hierarchies are considered in the 
adiabatic oscillation limit. 
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Figure 16. Radial profiles of the neutrino emissivities for the pair annihilation and the plasmon decay in the 15 Mq just prior to collapse. 






